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Abstract. We describe the embedded resolution of a quasi-ordinary sur- 
face singularity (V, p) which results from applying the canonical resolution 
of Bierstone-Milman to (V,p). We show that this process depends solely on 
the characteristic pairs of (V, p), as predicted by Lipman. We describe the 
process explicitly enough that a resolution graph for / could in principle be 
obtained by computer using only the characteristic pairs. 



INTRODUCTION. The Jungian approach to resolving the singularities of an 
embedded surface V C C 3 begins with a projection it : V — ► C 2 with (reduced) 
discriminant locus A C C 2 . Let a : (M, ct~ 1 (A)) — ► (C 2 , A) be an embedded 
resolution of the plane curve A. Thus, M is smooth (of dimension 2) and er _1 (A) 
has only normal crossings. The bimeromorphic map a and the projection n induce 
a bimeromorphic map a' : V — V x C 2 M — > V, and a projection tt' : V' — ► M 
whose discriminant locus has only normal crossings. We have the following diagram. 



V 



M 



V 



A 



pro- 
laryfl 



This leads us to consider surface singularities in C 3 admitting a local finite pro- 
jection to C 2 with normal crossings. Such singularities are called quasi-ordinary. 
Lipman and others have studied this kind of singularity in detail (c/. fLU , L2 
and @). Each quasi-ordinary singularity is equipped with a set of pairs {(A, //)} of 
non-negative rational numbers, called characteristic pairs. These pairs determine 
much of the geometry and topology of the singularity. 

One approach to understanding a quasi-ordinary singularity is to study its reso- 
lution. In [LI |, Lipman described a definite procedure for resolving a quasi-ordinary 
surface singularity (V,p) embedded in C 3 (for an outline of this procedure, see [L2|, 
pages 164-165 and pages 169-171). He shows that at each stage of the resolu- 
tion, the singular locus remains quasi-ordinary, and the characteristic pairs at each 
stage are determined by those of the previous stage and the process employed. In 
this paper we study the embedded resolution of a quasi-ordinary surface singular- 
ity. Recently, several constructive algorithms for canonical embedded resolution of 
singularities have emerged; see Bierstone and Milman ([BM1], [BM2|), Villamayor 



1 Higher dimensional quasi-ordinary singularities are similarly defined. In this paper we only 
consider surface singularities in C 3 . 
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( | EV | ) , and Moh (Q)- In this note we apply Bierstone and Milman's algorithm. 



Specifically, we replace Lipman's procedure by applying Bierstone and Milman's 
algorithm for canonical embedded resolution, and show that at each stage of the 
algorithm, the center to blow up is determined by the characteristic pairs and the 
process employed. Thus we establish a canonical embedded resolution for quasi- 
ordinary surfaces; the fact that the data of the resolution depend only on character- 
istic pairs will be used in a sequel to prove the existence of simultaneous embedded 
resolution for equisingular families. 

The authors are thankful for several conversations and/or correspondences with 
E. Bierstone, G. Kennedy, J. Lipman, O. Villamayor and S. Encinas. A partial 
version of these results were obtained earlier in terms of the algorithm of jv|, 
|EV[ . We ultimately found it easier to use the more explicit bottom- up approach 
of Bierstone and Milman's treatment of canonical resolution, but we learned this 
subject equally from E. Bierstone, O. Villamayor and S. Encinas. We gratefully 
acknowledge our debt to them. We are told that S. S. Abhyankar has presented a 
proof of our main result in a new appendix to the second edition of his book 0, 
which we have not seen, but which may not contain an explicit recipe in terms of 
the characteristic pairs. 

1. Preliminaries 

1.1. Quasi-ordinary Singularities. A quasi-ordinary singularity is an analytic 
germ (V, p) of dimension d which admits a finite map of analytic germs 7r: (V,p) — > 
(C d ,0) whose discriminant locus A (the hypersurface in C over which 7r ramifies) 
has only normal crossings. We consider a quasi-ordinary surface singularity embed- 
ded in (C 3 , 0) here. In this case, we can choose local coordinates x, y, z such that 
Tv(x,y,z) — (x,y) and such that (V,p) is defined by an equation 

f(x,y,z) = z m + c m - 1 [x,y)z m ~ 1 + ■ ■ ■ + co[x,y) = 

where the Ci are power series. Then {V,p) being quasi-ordinary means that the 
discriminant D(x, y) of / (considered as a polynomial in z) has the form 

D(x,y) = x a y b e(x,y), e(0,0)^0. 

It is known that the roots of / are fractional power series. Let £ = H(x 1 ^ n ,y 1 ^ n ) 
be one root of / (H(x, y) a power series). The other roots (the conjugates of Q can 
be expressed as 

= Hiix 1 '" *,y 1/n ) = ff( Wil z 1/n ,cW /n ) 

with the u>ij the n-th roots of unity. By letting d = £, the polynomial / can be 
written as 

f(x,y,z) = z m + c m -i(x, y)z m ~ 1 -\ \-co(x,y) 

m 

j=i 

so that the Ci are the elementary symmetric functions of Q . We call £ a parametri- 
zation of (V,p). 

Since x a y b e(x, y) = D(x, y) = ]X>y (Ci - G)> we have 

- o = * u/ V /n ^(* 1/n ,y 1/n ), e«(o,o) + o 
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for some u and v (depending on i, j). The fractional monomials My = x u / n y v ' n so 
obtained are called the characteristic monomials of C, and the exponents (u/n, v/ri) 
are called the characteristic pairs of £. These pairs satisfy certain conditions (cf. 
|L2| , Proposition 1.5); for example, they are totally ordered by (Ai,/Xi) < (\2,p,%) 
if and only if x Al y /il divides x X2 y^ 2 (i.e. Ai < A2 and fix < 112)- It turns out that 



these pairs determine quite a lot of the geometry and topology of (V,p) (cf. [Ll|, 
L|, (O), and |§ 



Although a quasi-ordinary singularity may have different characteristic pairs 
which depend on the choice of Cj the characteristic pairs of a normalized para- 
metrization determine and are determined by the local topological type of (V,p) 
(cf. [Q]); that is, there is a set of characteristic pairs which are independent of the 
choice of £. 

Definition 1.1. A parametrization 

Q = x a ' n y h / n H{x 1 / n ,y l ' n ), 
where H (0, 0) 7^ is normalized if 

1. a and b are not both divisible by n, 

2. if a + b < n, then both a > and b > 0, 

3. labeling the characteristic pairs (\i, /J>i)i<i< s of C so that Ai < A2 < • ■ ■ < X s 
and Hi < (j, 2 < ■ ■ ■ < fj, s , we have (Ai, A 2 , . . . , A s ) > (pufJ®, ■ . . ,fJ- s ) (lexico- 
graphically) . 

Lemma 1.2. Any parametrization £ of a quasi- ordinary singularity (V,p) can be 
changed to a normalized one by the following type of change of coordinates: 

x' = x x 1 = z x' = y 

(i) y' = y (ii) y' = y (Hi) y' = x . 



p{x,y) 



z 



Proof. Let p(x,y) be the sum of the terms x a yP in £ with a, (3 both integers. By 
a change of type (i), we get a parametrization £' = £ — p(x, y). By the property of 
characteristic pairs (cf. p2| |, Proposition 1.5), the smallest characteristic monomial 
divides every term in So, C = x Xl y^ 1 H(x 1 ^ n , y 1 ^ 71 ) for some H. Therefore, we 
may assume that ( = x Xl y' 11 H(x 1 ^ n , y 1 ^), and 1 in ( |l.l| ) is satisfied. If 2 in (IT) is 



not satisfied, from the proof of the inversion formula (cf. |g| , Append ix) a change 



of type (ii) can take care of this. Upon a change of type (iii), 3 in (1.1) can be 



obtained. □ 



1.2. Resolution of Quasi-ordinary Singularities. In jLl| J. Lipman described 
a definite procedure to resolve a quasi-ordinary singularity (for an outline of this 
procedure, see [L2], pages 164-165 and pages 169-171). This procedure determines 
a tree of special transforms of a quasi-ordinary singularity (V,p). Each special 
transform in this tree is again a quasi-ordinary singularity. The center to blow up 
at each stage is determined by the characteristic pairs, and the characteristic pairs 
of the resulting special transform are determined by the pairs in the previous stage 
and the process employed (cf. [L2], Theorem 3.3). Let 

(2) ( = x^H(x 1 / n ,y 1 / n ) 1 H (0,0)^0 
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be a normalized parametrization of a quasi-ordinary singularity (V,p), where p = 
(0, 0, 0) is the origin. Roughly speaking, we blow up a permissible curve (an equi- 
multiple curve) whenever it is possible, and a point otherwise. The possible per- 
missible curve is the plane section cut out by coordinate plane x = (if A > 1) or 
y = (if /i > 1). By Bierstone-Milman's algorithm, however, we will at sometime 
blow up a point even though there is a permissible curve. 

The following lemma by Lipman shows how the characteristic pairs change af- 
ter blowing up the origin (quadratic transform) or a permissible curve (monoidal 
transform), (cf. |E^j, Theorem 3.3.) 

Lemma 1.3. (omit i = 1 if the corresponding pair consists of integers.) 



Transformation 








Pairs of resulting branch 


MONOIDAL TRANSFORMATION 










Center (x, z) 










Xi - 


- l,Mi 


Center (y, z) 










Xi, Pi — 1 


QUADRATIC TRANSFORMATION 










"Transversal Case 


" (Xi - 












Direction (1:0 


:0) 








Xi + fi 


i - 1, Mi 


Direction (0 : 1 


:0) 








Xi, Xi 


+ Mi - 1 


"Non- Transversal 


Case" 


(Ai +/xi < 1) 










Direction (1:0 


:0) 




Xi + 


(1 + 


Mi)(l- 


Ai) „ 1 + Mi j 






(1 + 


Mi 


Mi 


Direction (0 : 1 


:0) 




Mi + 


W- 


Mi) 1 + Aj x 
Ai 


Direction (0 : 


:1) 




Xi(l 


- Mi 


) + MjAi 


AiMi +/Xj(l — Ai) 






1 


-Ai 


- Mi 


1 - Ai - mi 



In the case of quadratic transformation, we will call the directions (1:0:0), 
(0:1: 0), and (0:0:1) the x-, y-, and z-charts respectively. 



2. Embedded Resolution 

Let f{x,y,z) define a quasi-ordinary singularity germ at the origin, and let 
{(Xi, fii)} be the collection of characteristic pairs for /. We describe the embedded 
resolution of f which results from applying the canonical resolution algorithm of 
to /. We show that this process depends solely on the characteristic pairs 
of /. We also describe the process explicitly; from our description the character- 
istic data could be fed to a computer to output a resolution graph. To describe 
the process we divide the problem into two natural stages. First we show how to 
reduce the multiplicity of / for the first time; we call this "the base case." Then we 
show how to reduce the multiplicity starting from a "general configuration." For 
this we define a configuration of exceptional divisors arising after the multiplicity 
has decreased, and show that it is stable. 

2.1. Outline of Bierstone-Milman's canonical resolution. The centers for 
blowing-up an arbitrary germ / are chosen by a process that can be briefly outlined 
as follows. Using the Weirstrass Preparation Theorem, / is written as a pseudo- 
polynomial of degree m in a variable z, where m is the multiplicity of /. It can also 
be arranged that the coefficient of the z m_1 term vanishes. Then the surface N\ = 
{z = 0} defines a "maximal contact-space" for /. The coefficients of / are regarded 
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as functions on N\, one dimension lower than the ambient dimension. From this 
collection, a (weighted) "higher order multiplicity" is defined and a function is 
selected which achieves this multiplicity, which sets up an induction on decreasing 
dimension. As the procedure is iterated, the accumulated exceptional divisors are 
factored out or included in the coefficient set, depending on when they occurred in 
the history of the resolution. Because the codimension cannot increase indefinitely, 
the process must end. It does so in one of two ways: 

1. At some stage, there are no more coefficient functions. Then the last max- 
imal contact space is contained in the vanishing locus of the previous set of 
coefficient functions. The final multiplicity is infinity, the order of vanishing 
of 0. Since the contact space is automatically smooth, it becomes the natural 
choice of center for the next blow-up. 

2. A multiplicity of zero is obtained. This occurs when the set of coefficient 
functions is generated by a product of exceptional divisors (which are fac- 
tored out before the multiplicity is computed; zero is the order of vanishing 
of 1). In the language of 0], jHj, the coefficient ideal is principalized. Each 
component of the maximum locus of the invariant is an intersection of ex- 
ceptional divisors, which is smooth by the inductive requirement of normal 
crossings. There is a rule, based on the order in history of the exceptional 
divisors, which determines a component to select as the next center. 

The details of this procedure are organized into an invariant inv f which consists 
of a finite sequence of non- negative rational numbers. The terms of the sequence are 
just the weighted multiplicities Vi described above and the numbers Sj of exceptional 
planes occurring in stages of history between which the invariant decreased. The 
numbers Sj reflect the "age" of exceptional planes; for example, s\ counts the oldest 
planes, those which occurred before the first multiplicity v\ decreased. 

Ordered lexicographically, inv^ is an upper-semi-continuous function on the 
points of {/ = 0}. Its maximum locus defines smooth components which can be 
ordered to defi ne can o nical centers for blowing up. All of these facts are explained 
and proved in JgMjJ , [ |BM2| . 

For quasi-ordinary singularities, we define a "general configuration of exceptional 
divisors" , and show that it is stable under every transformation that is needed in 
the resolution. (The transformations we use are blow-ups of points and curves, 
permuting coordinates, and collecting integral powers, i.e. the transformations used 
in Lemmas |l.2j and |l.3| ) . We also show that whenever the multiplicity decreases and 
the strict transform is normalized, the resulting collection of exceptional divisors 
constitutes a general configuration. 

The resolution will be presented in two stages; the first stage is the base where 
E, the collection of exceptional hyperplanes, is empty. The second stage is the 
case where the multiplicity has just dropped, and E is the general configuration. 
The first stage sets up the second; we will also use calculations in the first stage to 
introduce the machinery of Bierstone-Milman's canonical embedded resolution. 



2.2. Stage One. Let / = z m + - ■ -+cq(x. y) be a quasi-ordinary germ in normalized 
form, and let (A, /j,) be the first characteristic pair. Assume this is the beginning of 
the resolution and E = 0. By a series of canonically-chosen centers, we blow-up / 
until its multiplicity has been lowered. We distinguish two cases: (i) A + // < 1 (the 
"non-transverse case" according to Lipman) and (ii) A + fi > 1 (the "transverse 
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case"). It is easy to see that the multiplicity v\ of / is to (A + //) in case (i) and m 
in case (ii). 

Case (i) A + p, < 1 

Since / is normalized, A/i ^ 0. It follows immediately that the multiplicity of / 
at the origin is strictly greater than at any nearby point. Since V\ is the first entry 
in inv/, the invariant is maximized at the origin. Thus in this case, the canonical 
center is immediately determined by the first characteristic pair: the center is the 
origin. In particular it's not necessary to compute the entire sequence inv/ to 
determine the center in this case. 

We blow up the origin and examine three coordinate charts where the strict trans- 
forms are determined by x~ Ul f(x, xy, xz), y~ Vl f(xy, y, yz), and z~ Ul f(xz, yz, z), re- 
spectively. We will call them x-, y- and z-charts. (Lipman used the term (1:0: 0), 



(0:1:0) and (0:0:1) directions, respectively in |L2 .) 
In the x- chart we have 

x- Vl f(x,xy,xz) = x-" /1 {x m z m + --- + c (C l (x,xy)) 
= x m - Ul z m + --- + y Uim - unit. 



Using Lipman's formula (Lemma 1.3), the multiplicity has been lowered to v[ = pm. 
To restore / to (perhaps non-normalized) quasi-ordinary form, exchange variables 
y and z and apply the preparation theorem. We obtain 



/' = ^ m + --- + c'o(C' J ) 
(AV)- A ~ A_/i 1 



where the new data is: 



v \ = pm < (A + p)m = v\. 

Note also that p! > 1 and so /' is in the transverse case. /' is again a pseudo 
polynomial defining a quasi-ordinary singularity, but it may not be normalized: a 
branch may contain integer powers. To normalize /', replace z by z — p(x, y) where 
p is the (convergent) power series of all terms in £' with integral exponents. Assume 
for now that A', p! are not both integral, i.e. that /' does not lose a characteristic 
pair. It follows from the characterization of quasi-ordinary parametrization that 
x >- yP- \p(x,y). The final configuration is 

(3) /' = z m ' +--- + {x x 'y»') m -u, E = {(y)} 

and the multiplicity has dropped. 

The y-chart is exactly parallel. In the z-chart, several blow-ups of the origin 



are required (their number depends only on (A,/i)). Using Lemma 1.3 again, to 
and (A, /i) transform to m' = to — V\ and (A',/i') = , ); so long as 

the new (A',//) satisfies A' + fx' < 1, the multiplicity remains unchanged. The 
strict transform is automatically in quasi-ordinary form and the multiplicity is still 
maximized at the origin, which must be the next center. Iterating this process, we 
continue to blow up origins in the new z-charts until A'+/i' > 1. This last condition 

occurs after k = -j—^ — blow-ups, where [or] is the greatest integer function. The 

final configuration is 

/' = z m ' + ■■■ + x Xm y^ m ■ unit, E = {(z)} 



CANONICAL RESOLUTION OF A QUASI-ORDINARY SURFACE SINGULARITY 



7 



and m' = m(l — fc(A + //)) = v' \. By the definition of k, v'i < v\. Using Lipman's 
formula, 

(\> ,A = f^H Et\ = ( A V \ 

{ ^' \m"m'J yi-kiX + n)' l-k(X + fi)J ' 

So A' + n' > 1, and /' is in the transverse case. 

In x- or y-charts occurring after the first z-chart, the multiplicity drops after 
just one blow-up. The final configuration is the same as (||) except E also contains 
(z). 

Case (ii) A + [i > 1 

We have in this case v\ — m. We assume that A/i ^ 0. (The case where A/i = 
is simpler.) We assume that / is normalized and £ contains no terms with both its 
exponents integral. It follows that the coefficient of z rn ~ 1 in / is zero. Therefore in 
this normalized quasi-ordinary germ has a smooth maximal contact space 

given by the hyperplane defined by z = 0. We will illustrate the calculation of invy 
and the determination of centers for a few steps, and then summarize the resolution 
for this case of stage one. 

Proceeding inductively on dimension, let Hi be the collection of coefficients of 
/, regarded as functions on Ni = {z = 0}. More precisely, Hi is the set of pairs 
(c q ,m — q) where c q is the coefficient of z q . The number m — q is the weighted 



multiplicity assigned to c q . Following [BM1], define 



li 2 ■= mm <^ | (h, Hh) G ri\ 

where ^(h) is the order of h at O G iVx. While E — we define v 2 — ^2', v 2 is the 
"second multiplicity" at O. Now let 

Q-2 = {(hi,V2 ■ fih) I {h,Hh) G Hi}. 

Then there must be a pair (<?*,Ai g „) G Q2 such that = \i g , ■ From here we 

repeat the construction: we find a (linear) change of coordinates which leaves z 
unchanged such that a second variable defines a maximal contact space N2 C Ni 
for the function 5*. Restricting the coefficients of to N2, define H2 and 1/3. 



The salient points of this construction are: (see [BM1] and [BM2| for proofs and 
details) 

• The numbers i>i are analytic invariants of the singularity. 

• The construction at latter steps of the resolution is amended to reflect the 
presence of exceptional divisors. 

• The collection Qi can be replaced by an equivalent collection to make cal- 
culations simpler; an "equivalent" collection defines the same value for the 
invariant in a way that is stable after blowing-up. If (g, ji g ) G Q 2 satisfies 
n(g) = fig, and g splits into factors g — Yig™ 1 , then (g,fi g ) can be replaced 
by the pairs (gi,fJ- gi ), = fi>(9i), to obtain an equivalent collection. 

Now we complete the calculation of inv/ in the current setting of A + [i > 1 and 
E = 0. We have 

m 

f = z m + c m _ 2 (x, y)z m - 2 + ■■■ + co(x, y) = J](z - Q. 
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Each c q is a symmetric function of the Q. We have 

m 

c (x, y) = H Ci(x, y) = x mX y m ^ ■ unit, 
i=i 

and 

c q (x,y) = (x x yn m - 9 -e q (x,y), 

where e q is a fractional power scries, possibly a non-unit. This is a general phe- 
nomenon: for any quasi-ordinary singularity, cq "divides" the other coefficients in a 
weighted sense: hi "divides" hi if h^ 2 \ h% hl ■ So, c determines the next weighted 
multiplicity. Therefore, 

m(A + /i) 

"2 = M2 = = A + /i. 

m 

The pair (g,/J. g ) = {cq,v 2 ■ m) E Qi satisfies ^i(g) = [i g , and since c factors, 
we can replace it by the pairs (x, 1) and (y, 1). (Here we use the assumption 
A/i ^ 0.) The function (y, 1) already defines a second maximal contact space: 
N 2 = (y, z) C N\ = (z). As before, Hi is the collection of coefficient functions of 
the elements of Qi, written as polynomials in y. Thus (x, 1) G Hi, and all other 
functions (h, ^h) of TL 2 satisfy fj,(h) > [i^. Therefore, is 3 = fi 3 = 1. There are no 
more variables left to define coefficient functions, and ^4 is defined to be 00. We 
define Sj = for i = 1, 2, 3 and 

mvf = (^i,si;^2,s 2 ;^3,s 3 ;^4) = (m, 0; A + n, 0; 1, 0; 00). 

When invy ends in the value 00, the last contact space is chosen as the center. In 
this case, 7V 3 = (x, y, z) defines the origin. Note that the invariant depends only on 
the characteristic data of /, since m is determined by that data. 

The procedure for reducing the multiplicity falls into a handful of cases; we'll 
describe one case completely and then summarize the others. 

— Year One 

Blow-up the origin and examine the chart defined by x. The strict transform 
/' remains in quasi-ordinary form, assuming the multiplicity has not dropped. Let 
O' be the origin in this chart. Set E(0') = {H x } where H x = (x) is the excep- 
tional divisor. Define E 1 C E(0') to be the collection of exceptional hypcrplancs 
{H' I H' is the strict transform of H G E(O)} where O is the point that was blown 
up. Since E(0) = 0, i? 1 (0') = 0. In general, i? 1 (0') are those exceptional hyper- 
planes at O' which are pulled back from the earliest stage of resolution in which 
/ had the same multiplicity as it does now. E l+l (O r ) C E \ U' l q=1 E q is defined 
similarly, replacing the multiplicity by the first part of the invariant, namely the 
string hiVj_ 1 = (ui, s\; ■ ■ ■ The numbers Si are defined as the cardinalities of 
the sets E l . 

For clarity, we'll work through the machinery for calculating inv^ in the current 
chart; after this we give just the essential calculations. 

Let Q\ = {(f',m)} be the strict transform of /, which we've assumed retains 
the same multiplicity. Let T\ = G\ U where (E\ 1) = {(g,l) | g G E*} 

(Currently, E 1 = so Qi = From T\ we draw a function with /j,(g) = fi g , 

and use it to define N\. Of course, (/',m) is the only choice; since /' is already 
expressed as a normalized quasi-ordinary singularity, we may take Ni — (z). As 
before, Hi = {(c/, m — i)} where the c/ are the coefficients of /' as a polynomial in 
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■m(X' + u!) 

We see that \ii = = A' + //, where (A',//) is the first characteristic 



pair for /', which can be obtained from Lipman's formula (1.3) 



In general, [12 is modified to take exceptional divisors into account. Let (yn) = 
H 6 £1 = E \ E 1 and define piH = mm - where fiii(h) denotes the order to 

which yn factors out of h. Then v% := \ii — J]] A*2iJ- Define Z?2 := ][J yn^ H ■ 

He£i Hes t 
By construction, is the greatest common divisor of the elements of Hi that 
is a monomial in exceptional coordinates yn, H G £±. For each h 6 Hi, write 
h = D™ ■ g and /i g := fih ■ ^2- Then £2 is defined to be the collection {(g,H g )} for 
all (h,fih) S Hi as above, together with the pair (f™, (1 — v i) • m ) if ^2 < 1. This 
completes a cycle in the definition of inv f . 

Now calculate explicitly: <V = (x A y M ) m • unit is the principalizing element of 
Hi, and £\ = E \ E 1 = {H x } where H x — (x). Clearly /i// x = A', so — fx' and 
D2 = x x . We will assume that A > 1 and /i > 1; using Lipman's formula, \i = pi ^ 
so v<i > 1, and D2 is not included in C?2- We may exclude functions from Q2 which 
are divisible by some other function in C?2- We have Q2 = {(y M m , p'm)} ~ {(y, 1)}. 
The first stage of the invariant, invj^i, is (m, 0;/x') = (m, 0;/x). Before blowing up, 
invj^i = (m,0;A + /i), so invj^i has decreased. Therefore E 2 = {H x }, and £2 is 
defined to be £ x \ E 2 = 0. We have s 2 = 1, and «F 2 = £?2 U (£; 2 , 1) = {(x, 1), (y, 1)}. 
It's now clear that we may take (y,z) and (x,y,z) to successively define maximal 
contact spaces, and i/ 3 = 1, S3 = 0, = 00. We have inv/ = (to, 0; //, 1; 1, 0; 00). 
The last maximal contact space is the origin, which must be the next center to 
blow-up. 

— Year Two 

Blow up the origin, and again examine the x-chart. As before, the strict trans- 
form is a normalized quasi-ordinary polynomial. The assumption A > 1, ji > 1 im- 
plies the multiplicity cannot decrease by blowing up the origin, so v\ — m. It follows 
that E 1 = 0, and T\ = Gi contains only the strict transform, which we will again 
denote as /. Then as before Ni — (z), and Hi — {(co,to)} (we ignore all higher 
coefficient functions, which are divisible by Co). Again we have p% = A" where 
(A",/i") is the current first characteristic pair, and £1 — {H x } leads to fi2H x = A" 
and V2 = /i". Using Lipman's formula, \i" — (i and invj^i = (m,0;/i). This time 
invji does not drop, so E 2 — (the exceptional divisor from the previous stage is 
blown-away), and £ 2 = £1 \ E 2 = {H x }. Then T 2 = Q 2 = {{y^ m ,pm)} - {{y, 1)}. 
The next maximal contact space N 2 is defined by y restricted to Ni, i.e. iV 2 = (y, z). 
There are no coefficient functions contributed to H2 by the only function (y, 1) be- 
longing to Ti, so 1/3 = 00. We have inv/ = (to, 0; /x, 0; 00) at this stage. The next 
center is defined by N2, which is the x-axis. Note that by assuming fi > 1, the 
x-axis is equimultiple. 

— Year Three 

Blow-up the x-axis. In this case there is only one chart to consider, the y-chart. 
(The strict transform of / is smooth in the other chart.) The strict transform /' 
of / is y~ m /(x, y, yz). f is (up to multiplication by a suitable unit) a normalized 
quasi-ordinary polynomial, with first characteristic pair (A',//) = (A, fi — 1) where 
(A, /j,) was the corresponding pair in year two. To see the general pattern, we 
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continue to assume A', [n! > 1. Therefore the multiplicity m has remained constant. 
We see that E = {H x , H y } and E 1 = 0. Clearly T x = Gi = {(/',m)} ; and we 
take iVi = (z). The coefficients Tii are generated again by ((x A j/ M ) m ,m), but now 
£ 2 = E\E 1 = {H x , H y } implies D 2 = x x 'y»' '. Then [i 2 = A' + fi' and v 2 = 0. The 
invariant sequence terminates if = or oo. Thus invj = (m, 0; 0). 

When Vi = 0, the center for the next blow-up is constructed out of the exceptional 
locus. List the elements of E by the order in which each arose in history; this can 
be done globally, not just for the current chart. (Thus we consider E to contain n 
elements after n blow-ups, most of which play no role in the calculation of inv/ at a 
given point). Any subset I oi E can be ordered via the n-tuple (Si,5 2 . . . . S n ) where 
5i — if the ith element of E is not in /; n-tuples are ordered lexicographically. 

Returning to the determination of the center of blowing up, when v 2 = 0, define 
02 = {{D 2l 1)} and Sj nv = {peN 1 | ^ p (D 2 ) > 1}. 

In the current setting, 

Sinv = {peN 1 | V) > 1} = (H x n Nx) U (H y n iVi), 

the union of the x-axis and y-axis. The order of H x is (0, 1, 0) and the order of H y is 
(0, 0, 1). Thus H x precedes H y and the next center is C — fl H x = the j/-axis. 

The procedure now takes on a stable pattern. So long as A and \i both exceed 
1, we have inv/ = (m, 0;0). We are forced to continue blowing up axes, which 
lowers either A or \x by one each time. The corresponding divisor is blown away 
and replaced by a newer one, causing the choice of axes to alternate each time. 
Eventually one of A or fi is driven smaller than one, and the requirement n P (D 2 ) > 1 
automatically selects the other axis ever after, until both A < 1 and /j, < 1. If at 
this point we still have still have A + fi > 1, then the multiplicity still remains 
unchanged. But the condition (i p {D 2 ) = /z p (a; A y M ) > 1 now selects the origin as the 
center, and in fact the origin is the maximum locus for the multiplicity v\. We are 
forced to blow up the origin repeatedly until A + fi < 1. (Only x-charts and j/-charts 
need be considered.) The multiplicity drops precisely when this occurs. Finally, it 
is clear that all of these decisions are completely determined by the arithmetic of 
the original pair (A, p). 

Except for which charts (x or y) we choose to consider in the first two years 
and the assumptions on (A, /i), the preceding completely describes how to reduce 
the multiplicity starting from year zero in the transversal case. Other starting 
assumptions {e.g., A < 1, /i > 1) are entirely similar. 

2.3. Stage Two. Once the multiplicity drops, the algorithm begins again; all pre- 
vious history is forgotten. All exceptional divisors passing through the current point 
of interest are placed in E 1 (the set of "oldest" divisors). When the multiplicity 
drops for the first time, all the divisors are given by coordinate functions. In later 
stages this may not be true. The general phenomena occurring in later stages is 
completely illustrated by three examples. 

Example 2.1. Suppose / is normalized, X + p < 1, and E = E 1 = {H x ,H y }. This 
is the usual situation the first time the multiplicity is lowered when we begin in the 
transversal case. Also, note that A/x ^ since / is normalized. 

As before, A + fi < 1 immediately implies the origin has higher multiplicity than 
nearby points. We are required to blow up the origin. Examine the x-chart. As 
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before, the multiplicity drops immediately. The new first characteristic pair is 

and the strict transform of / has the form 

/' = z m x m - v + --- + y» m - unit. 

Note that fi' > 1 so we are again in the transversal case. We have E = E 1 = 
{H x , H y }. To normalize /', we first permute y and z, and, up to a unit factor, get 
the configuration 

(4) f = z» m + --- + x m -^y m -unit, E = {H X ,H Z }. 

But now /' may not be normalized; a quasi-ordinary branch may contain 
terms with integer exponents. Let p(x, y) be the power series consisting of these 



terms, and apply transformation (i) of Lemma 1.2. It follows that p(x, y) is divisible 

lz-p(x,y) 



by x x y^ . We now have H z _ p t xy \ 6 E 1 = E. Since N\ — {z — 0} can be chosen as 



the first maximal contact space, both (x x m y M m , to') and (p(x, y), 1) will belong 
to the coefficient set Hi, and by the divisibility relation between them, the latter 
function will have no effect on the calculation of inv^ . This relationship is preserved 
under blowing-up. Thus the divisor H z _ p ^ xy ) behaves as though p{x, y) = 0. 

Example 2.2. Suppose again that / is normalized and A + [i < 1 with E = E 1 = 
{H x , H y , H z _ p ( x ^} where x x y^ \ p(x,y). This case may occur when Example 2.1 
has lowered it's multiplicity a second time. 

Since / is in the non-transversal case, we must blow up the origin. In the z-chart, 
the divisor H z _ p ^ xy ^ is blown away and plays no further role. In the x-chart, the 
multiplicity drops immediately. As before, the strict transform /' requires two steps 
to be normalized: permute y and z, and then absorb integral powers in the quasi- 
ordinary branch. In the second step, we replace z by z—q(x, y) where x x y^ \ q(x, y) 

and (A',/i') = ( — , — ] . We need to see the effect of these changes on the 



elements of E . Of c ours e H x is unaffected, and H y becomes H z _ q ^ x y ^, which 



behaves as in Example 2T. On the other hand, H z - p t Xjy \ becomes Hy— p r XtZ —q(x,y))- 
By the divisibility properties of p and q, y^ | p(x,y) and x x y^ \ q(x,y). Since 
Ni = {z = 0} will again be the first contact space, it follows that the function 
y—p(x, z—q{x, y)) is equal to y up to a unit when restricted to Ni. Thus H y _ p ( x z _ q ) 
behaves like H y through out the next stage of the resolution. 



In Examples 2^ and 2/1, we implicitly assumed that Lipman's formula for (A', //) 



result in a new first characteristic pair which is not entirely integral. 

Example 2.3. Suppose / is normalized, A + /x<l, E = E 1 = {H x , H y , H z _ p ( xy )} 
with x x y^ | p(x, y), and (A, = (^, -M where 1 < a < to and a G N. 

As in Example |2.2| , blow up the origin and consider the direction defined by x. 
The multiplicity of /' drops, but this time (A', //) = (m — a — 1, m). (The new pair 
becomes integral.) When /' is normalized, the "first" characteristic pair (A', //) has 
to be absorbed in the transformation z — > z — q(x, y). Thus q(x, y) — x x y M • unit. 
The strict transform has one fewer characteristic pairs. (Eventually / has to lose 
all its characteristic pairs to become smooth.) 
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Just as before, H z becomes H y _ p ( xz _ q ( xy - ) - ) and the latter is generated by (y) 
once we restrict to Ni = {z = 0}. But while the transform H y — > H z _ q ( xy * ) as 
before, we now have /' = z m + • • • + (x X2 y^ 2 ) m ■ unit. The degree-zero coefficient 
of /', c ' = {x x ' 2 y^ 2 ) m ■ unit, no longer divides q(x,y). In fact the weighted 
"divisibility" relation is reversed: q(x,y) "divides" cq 1 . When we construct Hi, 
it is q that principalizes the functions and determines the invariant, until H z - q is 
blown away. Since the role of q in the process is determined by (A,/z), although 
(A', fj,') is integral and no longer counts as a characteristic pair, the process is still 
determined by the original characteristic data. We sketch the next phase of the 
resolution. The invariant for the configuration (f',E) is (m', 3; 1,0; 1,0; oo) and 
the next center is the origin. Blow-up the origin and consider the x direction; 
the multiplicity can't drop but H x is blown away and replaced by a new divisor. 
The invariant records this as (m 1 , 2; 1, 1; 1, 0; oo), and the origin is the center again. 
Blow up the origin and consider the y-direction; now H x and H y are "new" and 
belong to S\. The calculation now gives Hi — {(co,m'), (q, 1)} and D2 = x x y^ , 
where q(x,y) = x x y^ ■ unit. We have ^2 = 0, and inv/ = (m',l;0). D2 now 
determines the centers for blowing-up by S = {p | ^i p (_D 2 ) > 1} and the ordering 
of the divisors in Si = {H x ,H y }. So we are required to blow up the x- and y- 
axis alternately until A' < 1 and [i! < 1; since (A',/u') are integers, we continue 
until A' = fi' = and H z _ q i x>y -s is blown away. The relation q{x,y) \ c (x,y) is 
maintained though-out, so when H z _ q t x y \ is blown away we have Hi = {(co, m')}. 
Then since Co(x, y) = x X2 y^ 2 - unit, the second characteristic pair (A2, /U2) takes over 
the process. 

Definition 2.4. General Configuration. Let / be a normalized quasi-ordinary sin- 
gularity, and let (A, fi) be the first characteristic pair. Furthermore, suppose / was 
obtained by applying the Bicrstone-Milman algorithm to / with characteristic pairs 
A = {(A*, /i;)}. Let 

E = E 1 C {0 - q(x, y)), (y - q{x, z - p(x, y))), (x)} 

where 

y - q(x, z - p(x, y))\ z = y- unit, 

and x x y fJ - | q(x, y), or q(x, y) — x a y h -maX and x a y b | x x y^, and where (a, b) depends 
on A. Then E is called a general configuration. 

Theorem 2.5. At any stage in the resolution when the multiplicity of f(x, y, z) 
has dropped, and f is expressed as a normalized quasi- ordinary singularity, the 
exceptional divisors constitute a general configuration. 

Proof. The theorem is trivial in the base case. We have shown it to be true at 
the end of the base case, and in representative examples of the general case. A 
complete proof involves checking the calculations for each possible subset of the 
general configuration, and seeing that in each case the final set of divisors (once / 
is normalized) is again a general configuration. The calculations are similar in all 
cases to the three examples we presented above. □ 

Since the controlling function in every case depends only on characteristic data, 
we have 
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Theorem 2.6. The Bier stone- Milman canonical resolution algorithm for a quasi- 
ordinary singularity (V,p) depends only on the (normalized) characteristic pairs. 

Example 2.7. We will find locally an embedded resolution for (V,p) defined by 
f(x, y, z) = z 3 + x 2 y A . The function / has one characteristic pair (2/3, 4/3). 

— Year 

In the beginning, infj(O) = (3, 0; 2, 0; 1, 0; oo) and the center is the origin O. 
Blow up O and consider the x-chart. 

— Year 1 

f = z 3 + x 3 y 4 , inf/(0) = (3, 0; 4/3, 1; 1, 0; oo) and E = {H x }; the center is again 
the origin. Blow up the origin and consider the y-chart. 

— Year 2 

f = z 3 + x 3 y 4 and inf/(0) = (3, 0; 0) The centers are now determined by D% = 
xy 4 ! 3 and the ordering of the divisors E = {H x ,H y }. Since H x is older, we first 
blow up H x n N\ which is the y-axis. 

-Year 3 

f = z 3 + y 4 , inf/(0) = (3, 0;0), D 2 = y 4 ' 3 - Blow up the x-axis, and consider 
the y-chart. 

— Year 4 

f = z 3 +y. The multiplicity has dropped, and the strict transform is now smooth, 
but the divisors do not have normal crossing with /. If we "normalize" / as a quasi- 
ordinary singularity, we see that / has lost its only characteristic pair: the previous 
characteristic pair becomes (0, 1/3), which becomes integral when we normalize by 
the transformation y <-> z, z — > z + y 3 . We have f — z, and E = E 1 = {H x , H z _ y s}. 
Since the multiplicity just dropped, we begin again, regarding the current state as 
year 0. We have inf/(0) = (1, 2; 1, 0; 3, 0; oo) and the center to blow up is the origin 
O. Blow up and consider the x-chart. 

— Year 1 

In the x-chart we have f = z and E — {H x , H z - X 2 y 3 }- Since of course the 
multiplicity cannot drop, the new divisor H x does not belong to E . We have 
inf/(0) = (1, 1; 3, 1; 1, 0; oo) and the next center is the origin. Blow up and consider 
the y-chart. 

— Year 2 

f = z and E = {H x , H y , H z _ x i y i\. Now we have E 1 = {H z _ x 2 y i} and E\ = 
{H x , H y }, which leads to inf/(0) = (1, 1; 0) and D 2 = x 2 y A . All remaining centers 
are determined by D 2 , the exponents of which came from the integral pair (0,3) 
and the transformations of Lemma |l.2j occurring in the last two steps. As in the 
last cycle, we blow up the y-axis. 

-Year 3 

f = z and E 1 = {H z _ xy i}, and invj(O) = (1,1; 0). Blow up the x-axis and 
consider the y-chart. 

— Year 4 

f = z and E 1 = {H z _ xy 3}. Blow up the y-axis and consider the x-chart. 
-Year 5-8 
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In year 5, / = z and E 1 = {H z _ y z}. Now D2 = y 3 and only the x-axis is selected. 
Blow up the x-axis three more times. In year 8, / = z and E = {H x , H y , H z -i}. 
We have infy(O) = (1,0;0). The invariant inf/ is now locally constant everywhere 
on /, and we have achieved an embedded resolution. 
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